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Summary 


The research during the third reporting period has focused on fixed order robust control 
design for hypersonic vehicles. A new technique has been developed to synthesize fixed 
order Hoo controllers. A controller canonical form is imposed on the compensator struc- 
ture and a homotopy algorithm is employed to perform the controller design. Various 
reduced order controllers are designed for a simplified version of the hypersonic vehicle 
model used in our previous studies to demonstrate the capabilities of the code. However, 
further work is needed to investigate the issue of numerical ill-conditioning for large order 
systems and to make the numerical approach more reliable. 
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Introduction 


The objective of this research is to address the issues associated with the design of ro- 
bust integrated flight control systems for future hypersonic vehicles with airbreathing 
propulsion systems. It is anticipated that such vehicles will exhibit significant interac- 
tions between rigid body (airframe) dynamics, structural dynamics and engine dynamics. 
The uncertainty in the initial dynamic models developed for these vehicles will also be 
high. The main reason that highly interactive uncertain dynamics are to be expected is 
that scramjet engines will be the primary source of propulsion at hypervelocity speeds. 
Wind-tunnel testing as a result will be limited, and it will be necessary to gain experience 
in actual flight testing of such vehicles. This means that initial flight control system design 
efforts will rely more heavily on theoretical and computer based models, than has been 
the case for subsonic and supersonic aircraft. Also, propulsion system sensitivity to angle 
of attack variations and structural vibrations will lead to highly interactive dynamics. 


In this study, the current major research issues from a flight control viewpoint are: (1) the 
development of models that are representative of the interactive dynamics that can occur 
in such vehicles, (2) the development of representative uncertainty models for these dy- 
namics and (3) the development of practical approaches to designing multivariable flight 
control systems that guarantee adequate performance in the presence of uncertainty. The 
research done during the third reporting period has been focusing on item (3). 

The hypersonic vehicle model used in this study [1] neglects both the effects of angle 
of attack variations on propulsion system performance and of elastic body bending. In 
the framework of robust control design these effects are treated as uncertainties. This 
study has resulted in uncertainty models capturing the individual characteristics of these 
phenomena. The models were developed and incorporated into a control design structure 
which evolved from an earlier study described in Ref. [2]. A variety of robust controllers 
were designed utilizing Hqo and //-synthesis techniques and the sensitivity of achievable 
robust performance to the introduced uncertainty levels were investigated. A thorough 
description as well as a comprehensive discussion of the results is given in Ref. [3]. While 
demonstrating the feasibility of modern control techniques for hypersonic control system 
design, the results in Ref. [3] also indicate that a significant drawback of these meth- 
ods will be the large order of the resulting controller. When implemented, large order 
controllers can create time delays which may be undesirable. Furthermore, it was shown 
that model order reduction techniques to reduce controller order do not guarantee robust 
performance. 


These results demonstrate the need for fixed order controller design. A technique to 
synthesize H <*, controllers with a constraint on controller order is developed in Ref. [4] 
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which was submitted to the 1994 AIAA Guidance, Navigation, and Control Conference 
in Scottsdale, AZ, (see Appendix). Therefore, this report provides only a brief summary 
of the research performed and focuses on a discussion of future research. 


2 Fixed Order Controller Design 

The objective of fixed order design is to constrain the order of the controller a priori in 
the design process and to synthesize a robust controller subject to this constraint. In our 
case, a controller canonical form was imposed on the compensator structure. This leads 
to a closed loop formulation of the problem with a minimal number of parameters and 
converts it to a static gain output feedback problem. 


In Ref. [5] a conjugate gradient method was used for the numerical solution of this prob- 
lem. Besides the issue of slow convergence near the optimum, this method requires an 
initial starting guess which is sometimes hard to find. Thus, in our approach a homotopy 
algorithm is used to simplify this procedure. A homotopy starts with a known or easily 
calculated solution to a simple problem and deforms it gradually to the solution of the 
desired, but more complex problem while pursuing a so-called homotopy path. Here, the 
solution to the simple problem is a reduced order H 2 controller designed for low control 
authority. Using a two-step homotopy procedure, the desired fixed order H ^ controller 
for the original formulation is obtained. A detailed derivation of the problem formulation 
as well as an introduction into the homotopy algorithm used is given in Ref. [4]. 


The developed algorithm to synthesize fixed order H <*, controllers was applied to the 
Winged- Cone Configuration used in our earlier studies [1]. However, in order to perform 
initial calculations a simplified version of the ’’full-scale” model developed in Ref. [3] 
was used. Hence, the achieved results demonstrate the capabilities of the developed algo- 
rithm but do not necessarily represent realistic responses. Our future work will eliminate 
this shortcoming. Also, the current version of the algorithm is sensitive to numerical 
ill-conditioning for large order plants such as the one used in this study. Currently, re- 
search is under way to robustify the algorithm and to make the numerical approach more 
dependable. 


The uncertainty models developed in Ref. [3] representing propulsion system perturba- 
tions, flexible body bending, and uncertainty in control effectiveness are not included in 
these early stages of fixed order controller design. Since an design does not take into 
account the structure in the uncertainty, such a design would be too conservative in order 
to achieve robust performance. This is illustrated in Fig. 1 where the full order /z design 
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from Ref. [3] is compared to an equivalent full order Hoo design. The uncertainty levels 
are 25% uncertainty in CMa due to propulsive effects, uncertainty in the angle of attack 
caused by elastic deformation with the first body bending mode at 18.5 rad/sec, and 10% 
uncertainty in control effectiveness in both control channels. Clearly, the H <» design does 
not achieve robust performance while the // design does. This demonstrates the necessity 
to develop a fixed order robust control design technique that considers the structure in 
the uncertainty. 



Figure T. Robust performance bounds: design vs. fi design. (Also given: NP bounds 

(dash) and RS bounds (dash-dot) for /i design) 


3 Future Research 

The results achieved with the methodology described in the previous section are quite en- 
couraging and the approach has the potential for further development. The need for fixed 
order robust control design for structured uncertainty is obvious, but no such technique 
presently exists. Our future work will focus on the development of such a methodology. 


One possibility to pursue this task would be to replace the full order H c 0 design step 
in the ^/-synthesis procedure with the fixed order technique. Since //<*, controller design 
is a subproblem when designing for robust performance with structured uncertainty, the 
fixed order technique has the potential to constrain the order of the controller which is 
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normally subject to significant increases in the //-synthesis procedure. This “brute force” 
technique would adopt the D-K iteration procedure including the rather awkward step 
of curve fitting to find a rational transfer function representation of the D-scales . 


Currently, research is under way to circumvent this curve fitting phase in order to generate 
a completely automated design procedure for robust controller synthesis. In Ref. [6] the 
D-scales are selected to be constant instead of frequency dependent. This eliminates the 
need for curve fitting, and, moreover, allows performing a simultaneous optimization of 
D-scales and controller making a D-K iteration superfluous. Naturally, constant D-scales 
are more conservative than frequency dependent ones. Dynamic D-scales of a prespeci- 
fied order are used by the same authors in Ref. [7]. A state space representation of the 
dynamic D-scales is formulated and the optimal D-scales are computed. An attempt to 
simultaneously solve for optimal dynamic D-scales and an optimal fixed order controller 
using the formulation in Refs. [6] and [7] was not undertaken and would have lead to 
a large system of coupled equations. The numerical solution of such a system is rather 
difficult or may even be impossible. The canonical representation with its minimal num- 
ber of parameters has the potential to greatly simplify the formulation of this problem. 
Furthermore, the optimization procedure used in Refs. [6] and [7] would naturally fit into 
the general homotopy framework used in the 7/oc design. This issue will be investigated 
further. 


A different approach to bypass the curve fitting step is pursued in Refs. [8] and [9] where 
the authors use Km-synthesis for robust controller design. Employing a bilinear sector 
transform, the D-scales can be replaced by generalized Popov multiplier matrices M of 
a prespecified order. Similar to the D-K iteration in //-synthesis , a M-F iteration is 
performed where alternately full order Hoo controllers ( F ) and optimal multipliers (M) 
are computed. Embedded in this procedure is a fixed order multiplier optimization that 
replaces the curve fitting. The entire optimization is not simultaneous as in Ref. [6] and 
does not yield reduced order compensators, but it allows for controller synthesis using 
dynamic multipliers. Another major advantage of this technique is that both complex 
and real parameter uncertainty can be treated. The capacity to model variations in real 
parameters is considered a primary contribution that reduces conservatism in the design 
even further. Improvements in this area are currently the subject of expansive research 
activities [10] - [13]. In A r m -synthesis real uncertainty representation is accomplished by 
the utilization of the generalized Popov multiplier which allows for the introduction of 
phase information into the design. This technique is comparable to the approaches taken 
in Refs. [12] and [13]. It is not yet clear how the generalized Popov multiplier formula- 
tion can be included in the desired homotopy /fixed order design framework. However, it 
offers an interesting perspective especially under the aspect of including real parameter 
uncertainty. 
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An in depth application to the hypersonic vehicle model used in our previous work will 
be studied to evalute the developed techniques, and to investigate the critical modeling 
and design issues related to fixed order robust control design for this class of vehicles. 


References 

[1] J.D. Shaughnessy et al. “Hypersonic Vehicle Simulation Model: Winged-Cone Con- 
figuration”. NASA TM 102610, November 1990. 

[2] I.M. Gregory et al. “Hypersonic Vehicle Control Law Development Using Hoo and 
^-Synthesis”. In Proceedings of the 4 th AIAA International Aerospace Planes Con- 
ference , Orlando, FL, December 1992. 

[3] H. Buschek and A.J. Calise. “Robust Control of Hypersonic Vehicles Considering 
Propulsive and Aeroelastic Effects”. In Proceedings of the 1993 AIAA Guidance, 
Navigation and Control Conference , Monterey, CA, August 1993. 

[4] H. Buschek and A.J. Calise. “Fixed Order Robust- Control Design for Hypersonic 
Vehicles”. Submitted for presentation at the 1994 AIAA Guidance, Navigation, and 
Control Conference , Scottsdale, AZ, August 1994. 

[5] G.D. Sweriduk and A.J. Calise. “Robust Fixed Order Dynamic Compensation: A Dif- 
ferential Game Approach”. In Proceedings of the 1993 IEEE Conference on Aerospace 
Control Systems , West Lake Village, CA, May 1993. 

[6] W.M. Haddad, E.G. Collins, and R. Moser. “Structured Singular Value Controller 
Synthesis Using Constant D-Scales without D-K Iteration”. Preprint, Sept. 1993. 

[7] W.M. Haddad, E.G. Collins, and R. Moser. “Computation of the Complex Structured 
Singular Value Using Fixed Structure Dynamic D-Scales”. Preprint, Sept. 1993. 

[8] M.G. Safonov and R.Y. Chiang. “Real/Complex R m -Synthesis without Curve Fit- 
ting”. Control and Dynamic Systems , Vol. 56, pp. 303-324, 1993. 

[9] M.G. Safonov, J.H. Ly, and R.Y. Chiang. “^-Synthesis Robust Control: What’s 
Wrong and how to Fix It?”. In Proc. 1993 IEEE Conference on Aerospace Control 
Systems , West Lake Village, CA, May 1993. 

[10] M.K. Fan, A.L. Tits, and J.C. Doyle. “Robustness in the Presence of Mixed Para- 
metric Uncertainty and Unmodeled Dynamics”. IEEE Transactions on Automatic 
Control, Vol. 36, No. 1, pp. 25-38, January 1991. 


6 


[11] P.M. Young, M.P. Newlin, and J.C. Doyle. “/x Analysis with Real Parametric Uncer- 
tainty”. In Proceedings of the 30 th Conference on Decision and Control , Brighton, 
England, December 1991. 

[12] W.M. Haddad and D.S. Bernstein. “Parameter-Dependent Lyapunov Functions and 
the Popov Criterion in Robust Analysis and Synthesis”. Submitted to IEEE Trans- 
actions on Automatic Control. 

[13] J.P. How. Robust Control Design with Real Parameter Uncertainty using Absolute 
Stability Theory. Ph.D. thesis, Massachusetts Inst, of Technology, Jan. 1993. 


7 


Appendix 



Draft Paper for AIAA Guidance, Navigation and Control Conference 
Scottsdale, AZ, August 1-3, 1994 


Fixed Order Robust Control Design for 
Hypersonic Vehicles 

Harald Buschek* and Anthony J. Calise** 

Georgia Institute of Technology 
Aerospace Engineering 
Atlanta, GA 30332 


Abstract 

Control system design for hypersonic vehicles using modern control theory techniques 
generally leads to high order controllers. The time delays which are created when these 
controllers are implemented may not be acceptable. A technique to constrain controller 
dimension apriori in the design process employing a controller canonical realization and 
a homotopy algorithm is used to design fixed order controllers. Various reduced order 
controllers are designed for a hypersonic vehicle model accelerating through Mach 8. A 
comparison with a full order design is conducted and the results indicate that the devel- 
oped fixed order design technique is suitable for hypersonic applications. 
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1 Introduction 


At the dawn of the next millenium, hypersonic flight is one of the most demanding fron- 
tiers challenged by aerospace scientists. Designing a vehicle that has horizontal take-off 
and landing capability, is able to perform hypersonic atmospheric flight and reaches a typ- 
ical low earth orbit involves considering a variety of different effects ands requirements. A 
crucial component of such an aerospace vehicle is the control system which stabilizes the 
vehicle and ensures precise tracking of pilot commanded inputs. This task will be further 
complicated by a strong interaction between aerodynamics, structure, and the propulsion 
system. 

It has been shown in Ref. [1] that modern control theory using H^&ud //-synthesis 
techniques is well suited for addressing multiple uncertainty sources in hypersonic flight 
control design. However, a significant disadvantage of these techniques is that the result- 
ing controller is of the same order as the generalized plant. Frequency dependent weights 
have to be included in the design framework in order to achieve the desired performance 
and robustness characteristics. Thus, the order of the generalized plant is increased re- 
sulting in high order controllers. When implemented, large order controllers can create 
time delays which may be undesirable. One solution to this problem is to use model 
order reduction on the controller realization. This technique, though, does not consider 
the properties of the closed loop system when reducing the order of the controller and 
therefore robustness properties are not guaranteed. Another approach to this problem is 
to constrain the order of the controller apriori in the design process. 

In this paper, a controller canonical form is imposed on the compensator structure to 
constrain the controller dimension. Necessary conditions for an optimal /^controller are 
derived using a differential game formulation. A homotopy algorithm is used to continu- 
ously deform the solution of an H 2 problem formulation which serves as starting point to 
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the solution of the desired problem formulation, i.e. an optimal fixed order controller. 
This technique is applied on a hypersonic vehicle accelerating through Mach 8. Several 
fixed order Hoo controllers are obtained and their performance is compared to a full order 
design. 


2 Fixed Order Compensators for the iiooProblem 


For a standard control problem, the generalized plant is given by 


X = 

Ax + B\w + B 2 u 

(1) 

z = 

C\X + D 12 U 

(2) 

y = 

C 2 x + D 2 \W + D 2 2 u 

(3) 

where x £ 3?” is the state vector, 

w £ %l nw is the disturbance vector, 

u £ 9f i nu is the 


control vector, 2 £ 3J nz is the performance vector, and y £ $R ny is the observation vector. 
It is assumed that 

. (A,Bi,Ci) is stabilizable and detectable 

• (A,B 2 ,C 2 ) is stabilizable and detectable 

• Du has full column rank 

• D 2 i has full row rank. 

A general compensator for this system is . 

x c = A c x c + B c y 

XL — C cX c 


(4) 

(5) 
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where x c € 5R nc is the state vector of the controller the dimension of which can be speci- 
fied. Fig. 1 illustrates this design framework. 


The drawback with this general controller formulation is that the problem is over- 
parametrized. To avoid this problem of overparametrization, a canonical form description 
is used for the controller [2]. In controller canonical form the compensator is defined as 


P°x c + N°u c - N°y 

(6) 

-Px c 

(7) 

-Hx c 

(8) 


where x c <E 3R nc and u c € $ ny - P and H are free-parameter matrices, and P° and N° are 
fixed matrices of zeros and ones determined by the choice of controllability indices v, as 
follows: 


P° = block diag{i* (9) 

0 1 ••• 0 

i = (10) 

; i 

0 o 

N° = block diag{[0...01]f xi/i } i = 1, . . . , ny (11) 

The controllability indices must satisfy the following condition: 

5>i = nc i = l,...,ny (12) 

1 

Fig. 2 shows the structure of such a controller. Similarly, a compensator in observer 
canonical form can be constructed. However, in this paper only the controller canonical 
form is employed. 
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With this formulation, the compensator states can be absorbed into the generalized 


plant. Let 



The augmented system is defined by: 


x = 


A 0 

-n°c 2 p° 

Ax + B\w + B 2 u 


x + 


B 1 

-N°D 21 




C x 0 


X + 


D\2 0 


u 


y = 


C\X ■+■ D\ 2 u 
0 7 j x = C 2 x 

H 


u ~ — 


y = -Gy 


W + 


b 2 0 
-N°D 22 n° 


u 


(13) 


(14) 

(15) 

(16) 

(17) 


Eqs. 14-17 define a static gain output feedback problem where the compensator is repre- 
sented by a minimal number of free parameters in the design matrix, G. The augmented 
system is shown in Fig. 3. The closed loop system is given by 


x = (A — B 2 GC 2 )x + B\W 

= Ax + Bw (18) 

z = (Ci - D 12 GC 2 )x = Cx (19) 


In the /^problem, the objective is to minimize the oo-norm of the transfer function 
from disturbance inputs w to performance outputs z 


T zw = C{sl - A)~ l B . 


( 20 ) 


If z(s) = r zu ,(s)tw(s), then 

II^Hoo 


sup 

U , GZ/2 


l|Tzu>Hl2 

l|w||a 


sup \\T zw w \\ 2 . 
IMU<i 


( 21 ) 
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As described in Ref. [3], the optimization problem is to find 

mf{||r zu ,(G0||oo :GeG} := 7 * 


( 22 ) 


where G = {G G j R(" u + n ») xnc : A is stable}. A more practical design objective would be 
to find a G that insures ||T ztu ||oo < 7 for some 7 > 7*. Considering the performance index 

J 7 ( w, G) — E | J“(z T z ~ 7 2 w T w)dt} (23) 

where E{x( 0)} = 0 and £{z(0)x(0) r } = BB T , the H <*, problem can be formulated as a 
min-max problem 

J* = min max {J~(w, G)} (24) 

T G€G 7 ^6L 2 1 7V n 

where G 7 = {G € R( nu + n y) xnc : A is stable and H^Hoo <7}. In the min-max problem, 
the minimizing player acts first, and then the maximizer. Therfore, the loop is closed 
on u and w maximizes, whatever G might be. For any G G G 7 , a unique maximizing 
solution exists and is given by 

w — 'y~ 2 B T Q 00 x (25) 

where Qoo is the positive semidefinite solution of 

A t Q ao + Q 00 A + C T C + ^QooB^Qoo = 0 (26) 

Using the worst case disturbance in Eq. 25, the performance index can be rewritten as 

J 7 (G) = E { /°° x T (C T C - 7 - 2 Qoo5R T Qoo)xdf} . (27) 

Defining a distribution of initial conditions with zero mean and variance BB T , the objec- 
tive of minimizing H^Hoo using a fixed order controller can be formulated 

min {J 7 ((?) = tr{Q 0O 55 T }} (28) 

subject to the Riccati equation given in Eq. 26. In order to obtain the Ffoo-optimal 
compensator, the Lagrangian is defined as 

£(<?«,, L, G) = tr {QooBB 7 + ( A T Q 00 + Q^A + C T C + 'y- 2 Q 00 BB T Q 0 >)L} (29) 
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where L represents the Lagrangian multiplier. Matrix gradients are taken to determine 
the first order necessary conditions for an tfoo-optimal fixed order controller: 


<9L 

dQoo 

dL 

dL 

dL 

dG 


(A + + L(A + 7- 2 55 T goo) T + BB T = 0 (30) 

A t Q 00 + QooA + C T C + r'QooBBT'Q^ = 0 (31) 

Z{D T l2 DuGC 2 - DJA ~ B 2 Qoo)LCJ = 0 (32) 


Hence, three coupled equations have to be solved simultaneously to obtain the fixed order 
compensator which satisfies the constraint ||r 2U> ||oo < 7- 


Using this approach, fixed order i7oo-design can also be extended to fixed order /z- 
synthesis. Since Hoo controller design is a subproblem when designing for robust perfor- 
mance with structured uncertainty, the fixed order technique introduced above has the 
potential to constrain the order of the controller which is normally subject to significant 
increases in the /z-synthesis procedure. 


3 Controller Design Using a Homotopy Algorithm 

As demonstrated in the previous section, imposing a controller canonical form on the 
compensator structure provides a powerful tool for the design of fixed order controllers. 
Promising results have been obtained for the H 0 0 and the mixed problem in Ref. [3] where 
a conjugate gradient method was used. A disadvantage of this method is that convergence 
slows down near the optimum. Also, an initial starting guess for the compensator gain 
G has to be provided that stabilizes the closed loop system. In this paper, a homotopy 
method is employed to perform the fixed order controller design. A thorough discus- 
sion on homotopy methods and their utilization in fixed order controller design for H 2 , 
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H^as well as mixed ff 2 /#oo problems is given in Ref. [4]. Therefore, only a concise intro- 
duction into the principle of homotopy methods and their implementation in an -f/oc design 
procedure is provided in this paper. 

Homotopy (or continuation) methods, arising from algebraic and differential topology, 
embed a given problem in a parameterized family of problems. More specifically, consider 
sets U and Y € 3i n and a mapping F : U -*Y, where solutions of the problem 

F(u) = 0 (33) 

are desired with u€U and F(u)eY. The homotopy function is defined by the mapping 
H : Ux[ 0, 1] -*• such that 

H(u\, 1) = F(u) (34) 

and there exists a known (or easily calculated) solution, u 0 , such that 

H(u o ,0) = 0 (35) 

The homotopy function is a continuously differentiable function given by 

H(u(a), a) = 0, Va€[0,l] (36) 

Thus the homotopy begins with a simple problem with a known solution (35) which is 
deformed by continuously varying the parameter until the solution of the original problem 
(33) is obtained. 

For the i/oo problem two separate homotopy loops axe employed. An initial guess for 
the compensator gain Go is obtained by a full order, low authority H 2 design followed 
by order reduction of the compensator and transformation of the reduced order model 
into controller canonical form. This procedure was chosen because it has been observed 
numerically that order reduction techniques tend to work best for low authority LQG 
controllers [5]. A partial explanation of this phenomenon is given in Ref. [6]. The first 
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homotopy transforms the gain Go of the low control authority H 2 design to the gain 
G\ of a high authority H 2 design. This is done by gradually deforming the weights on 
control cost and measurement noise while proceeding along the the homotopy path until 
the desired problem formulation is recovered. A second homotopy is appended to per- 
form the i?oo design. The upper bound for ||T^||oo, 7, is reduced from some high initial 
value approximating the H 2 design to its minimal value for which a controller exists such 
that IIT^IU < 7. This procedure is illustrated in Fig. 4. In principle, it is possible to 
use a direct homotopy from G 0 to G while simultaneously deforming the weights and 7. 
However, the transition low-high authority always has to be completed in order to obtain 
the desired plant whereas the minimum value of 7 that can be achieved is usually not 
known beforehand. Thus, the proposed two-step procedure is more feasible for practical 
purposes. 


4 Performance Study for Hypersonic Vehicle 


The hypersonic vehicle model used in this study is the Winged-Cone Configuration de- 
scribed in Ref. [7]. In order to carry out control studies, a five state linear model of 
the longitudinal dynamics is used representing flight conditions for an accelerated flight 
through Mach 8 at approximately 86000 ft. The model has been normalized in the state 
variables velocity and altitude to avoid numerical ill-conditioning. State and control vari- 
ables are: 


vyioo 


velocity ((ft /sec)/ 100) 

a 


angle of attack (deg) 

<1 

— 

pitch rate (deg/sec) 

0 


pitch attitude (deg) 

h / 1000 


altitude (ft/1000) 
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u = 


(38) 


6e 


symmetric elevon (deg) 

Sq 


fuel equivalence ratio (-) 


The interconnection structure for controller design is shown in Fig. 5 and is based on 
the design in Ref. [1]. Inputs 1 to 5 introduce noise to the measurement signals velocity 
error K, altitude error h e , normal acceleration n z , pitch rate q, and pitch attitude 9 (out- 
put variables 7 to 11). Inputs 6 and 7 are command signals in V and h, inputs 8 and 9 are 
white noise signals introducing the turbulence components V g , a 3 , and a g . Inputs 10 and 
11 are the control signals in elevon and fuel equivalence ratio. Outputs 1 to 6 represent 
the performance variables. A description of the longitudinal model of the vehicle, the 
Dryden turbulence models and the actuator dynamics is given in Ref. [1]. Also given 
in Ref. [1] are uncertainty models introducing the effects of the propulsion system, the 
aeroelastic vehicle behavior and the uncertainty in control effectiveness. These effects, 
which have a considerable impact on performance, guidance and control characteristics of 
the vehicle, will be treated in the conference version of this manuscript. 


Fig. 6 shows a comparison of the maximum singular value plots of a full order design 
(12 controller states) vs. 9 th , 7 th , and 5 th order controllers. In the current formulation 
using the controller canonical form the minimal dimension of the compensator is 5 which is 
determined by the number of measurement outputs. The reduced order controllers come 
close, but do not fully recover the performance level of the full order controller which 
was obtained with the two-Riccati-equation solution described in Ref. [8]. This is due to 
the fact that the closer the homotopy approaches the optimal value of 7, the smaller the 
stepsize has to be in order to stay on the homotopy path. To reduce computation time, 
the code was stopped when the stepsize in 7 fell short of a prespecified tolerance, in this 
case A 7 mtn = 0.01. The values of 7 for which the design was stopped (jdesign) and the 
actual values of ||r 2tiJ ||oo obtained with this design (factual) are given in Table 1. 

For the 9^ order controller, the design value and actual value of 7 are very close. This 
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full order 

9 t/l order 

7 t,l order 

5 <h order 

7 design 

- 

0.4796 

0.8027 

0.7056 

7a dual 

0.3598 

0.4710 

0.7390 

0.6474 


Table 1: 7- values for different controller designs. 

indicates that the design is almost optimal although the level of the full order design is 
not achieved. For the 7 t/l and 5^ order controllers, the gaps between 7 design and 7 actual are 
larger implying that these designs can still be improved. In order to do so a finer stepsize 
would be required to enable the homotopy algorithm to track the solution closer to the 
optimal 7. An interesting feature is that the 5 th order controller achieves a lower 7 than 
the 7 tA order one. This can be caused by numerical difficulties the homotopy algorithm 
can encounter along its path forcing the stepsize to be reduced intermediately. This issue 
will be addressed in future research. 

An impression of how well the desired fixed order controllers perform is illustrated in 
Figs. 7 and 8. They show the velocity and altitude responses for simultaneous step inputs 
of V = 100 ft/sec and h = 1000 ft while encountering longitudinal and vertical atmo- 
spheric turbulence. It is evident that the responses for all reduced order controllers are 
almost as good and in most cases even better than the ones of the full order two-Riccati- 
equation solution. It is interesting to note that in Fig. 7 all systems with reduced order 
controllers exhibit a non-minimum phase behavior in the velocity response due to the 
pull-up maneuver to increase the altitude. This behavior may be eliminated by adjusting 
the weights applied to the error signals V e and h e . With the exception of the velocity 
response of the 5 state compensator all reduced order controllers tend to approach the 
desired steady-state values at least as good or better than the full order design. 
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5 Conclusions 


A technique for designing fixed order i/oo controllers has been developed and a control 
study of a hypersonic vehicle has been conducted. A controller canonical form descrip- 
tion was imposed on the compensator structure and a homotopy algorithm was used to 
design controllers of constrained dimension. A study of achievable nominal perfor- 
mance for various reduced order controllers demonstrates the capability of the developed 
algorithm and the feasibility of using fixed order control design for hypersonic aircraft. 
Reducing the allowable minimum stepsize in the homotopy will lead to even better results 
but will also slow down the procedure. Doing so as well as adjusting the performance 
weights is expected to reduce or the eliminate non-minimum phase characteristics ob- 
served in the velocity response. 


Important features typical for hypersonic atmospheric flight like propulsion system 
interactions and aeroelastic coupling effects are not yet considered. The final manuscript 
will concentrate on incorporating these features using the uncertainty models of Ref. [1]. 
Fixed order ^/-synthesis design will be used to address these multiple uncertainty sources. 
Additionally, the issue of further reduction of the controller dimension using an observer 
canonical form (in place of a controller canonical form) will be investigated, and issues 
related to controller architecture will be addressed. 
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Figure 4: Procedure for designing fixed order //^controllers. 
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Figure 5: Interconnection structure for controller design. 


16 








magnitude 



Figure 6: Maximum singular value plots for different controller designs. 



Figure 7: Velocity responses for different controller designs. 






Figure 8: Altitude responses for different controller designs. 
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